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Abstract 

The Huttner-Barnett model is extended to a magnetodielectric 
medium by adding a new matter field to this model. The eigenopera- 
tors for the coupled system are calculated and electromagnetic field is 
written in terms of these operators. The electric and magnetic suscep- 
tibility of the medium are explicitly derived and shown to satisfy the 
Kramers— Kronig relations. It is shown that the results obtained in 
this model are equivalent to the results obtained from the phenomeno- 
logical methods. . 
PACS number (s): 12.20.Ds 

1 Introduction 

Electromagnetic field (EM) is traditionally quantized by associating a quantum- 
mechanical harmonic oscillator with each mode of the radiation field in free 
space [H E]. This is achieved by introducing a general position and momen- 
tum and expressing the modes in terms of them. The transition from classical 
to quantum mechanical description takes place by treating the generalized 
position and momentum quantum mechanically and demanding equal-time 
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commutation relations (ETCR) between them. This scheme is called canon- 
ical quantization and is based on a Lagrangian density. 

In the presence of a linear polarizable medium canonical quantization is 
retained if the Hattner-Barnett model (HB) is used |3j. This model is based 
on the Hopfield model |4j and the matter is modeled by a harmonic field 
which represents the medium excitations. The inclusion of losses into the 
system can be done by inserting a reservoir consisting of a continuum of 
harmonic oscillators [5l [6] . Although this model is based on a microscopic 
model of matter, any microscopic parameter like coupling function between 
matter and EM field is not appeared in final results and in fact the results 
are written in terms of the electric susceptibility Xe which is a macroscopic 
parameter of the medium. 

Since the Hopfield model is written for a homogeneous polarizable matter 
and not for a magnetizable one, and since the HB model is based on this 
model, so we should extend the Hopfield model to a magnetizable medium. 
As mentioned, the final results of HB model do not depend on microscopic 
parameters, so we can change the microscopic relations such that we get the 
correct Maxwell equations in the magnetizable medium. 

Although by doing so, we may lose the physical interpretation of micro- 
scopic interaction, but at least we have a Lagrangian density from which 
correct Maxwell equations are obtained and a canonical quantization scheme 
is introduced. In addition, this model has a clear assumption about matter 
which can not be found in other macroscopic methods. 

The common method for the EM field quantization in the presence of 
matter is the macroscopic method (phenomenological method). The macro- 
scopic approach to EM field quantization begins with Maxwell equations and 
the loss is modeled by the Langevin force in the form of the noise current 
operators. A correlation function proportional to the imaginary part of the 
dielectric function is assumed for the later operations. In this way, a straight- 
forward calculation lead to the field expression in terms of the noise operators 
and the Green function of Maxwell equations [7] — [11]. In this method ETCR 
between the EM field operators are postulated and the verification of them 
justifies the validity of this approach. This model is recently extended to a 
magnetodielectric medium by adding a new noise operator for magnetization 
and assuming a correlation between them which is proportional to the imag- 
inary part of magnetic susceptibility [12] . In the conclusion of this paper, we 
will make a comparison between the results of our model and the macroscopic 
quantization model. 
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The paper is arranged as follows: In Sec. 2, the basic theory of electro- 
dynamics in a magnetodielectric medium is reviewed. In sec. 3, we extend 
the HB model to a magnetizable medium which is a special case of a mag- 
netodielectric medium, in fact, this section is an introduction to the more 
generalized case of a magnetodielectric medium. In sec. 4, we obtain a La- 
grangian density in a magnetodielectric medium from which we quantize the 
EM field. Finally, the main points are summarized in Sec. 5. 

2 Basic equations 

The classical Maxwell equations in a magnetodielectric medium are 

VxE(M) + ^ = 0. (1) 



VxH(r.t)-^ = 0. (2) 

The constitutive relations relate D and H to E and B. For isotropic linear 
medium, the relation between polarization and electric field is 

P(r, t) = £o f dt'xeir, t - t') E(r, t') dt', (3) 

J —oo 

where Sq is the permittivity of free space. The magnetization and the mag- 
netic induction field are related via a similar expression, 

M(r, t) = Ko f dt'xmir, t - t') B(r, t') dt', (4) 

J —oo 

where kq is the inverse permeability of free space. The electric and magnetic 
susceptibilities are defined in the frequency domain as 

X.(r,a;)=/ dr Xair,r) e'^\ (5) 
Jo 

where the index a could be e or m. The electric displacement D is expressed 
in terms of the electric field E and the polarization field P as 

D(r,i)=£oE(r,t) + P(r,i). (6) 
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The constitutive relation which defines the magnetic field H in terms of 
the magnetic induction B and the magnetization M is 



H(r,t) = ^oB(r,t)-M(r,t). (7) 

The real and imaginary parts of the electric and magnetic susceptibilities 
should satisfy Kramers-Kronig relations. This shows that the dissipative na- 
ture of medium is an immediate consequence of its dispersive character and 
vice versa. Therefore, in general we are dealing with a dissipative problem. 
For investigating a dissipative quantum system we should include a reser- 
voir in the process of quantization. After the quantization, the definition of 
polarization P, and magnetization M, in terms of the Fourier components 
should be reformed as 

POO 

P(r,t) = £o / M[Xe(r,a;)E(r,a;) + Pjv(r,^)]e^-* + //.c.}, (8) 
Jo 

and 

roo 

M{v,t) = Ko rfa;{[x™(r,^)B(r,a;) + Mjv(r,cu)]e^'^* + i/.c.}, (9) 
Jo 

respectively, where PAr(r,a;) and MAr(r,ci;) are polarization and magnetiza- 
tion noise operators and satisfy the fiuctuation-dissipation theorem 



3 Extension of the Huttner-Barnet model to 
a magnetizable medium 

Here, knowing the Lagrangian density of HB model, we change the form of the 
interaction term and find a Lagrangian density in a magnetizable medium. 
Let us write the Lagrangian as 

where 

= - |B^ (11) 

is the EM part which can be expressed in terms of vector potential A and 
scaler potential [/ (E = -A - VU and B = V x A). 

p up' 
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is the matter part, modeled by a harmonic oscillator field X of frequency cuq- 



Jo 




(13) 



is the reservoir part, consisting of a continuum of harmonic oscillators and 



is the interaction part which includes the interaction between the light and 
the magnetization field, with coupling constant a, and also the interaction 
between the magnetization field and the reservoir oscillator field with the 
frequency dependent coupling constant i/(a;). In this equation the interaction 
term between the EM and matter field is chosen such that it lead to correct 
Maxwell equations in a magnetizable medium. 

We make two assumptions about i>{uj): (i) the analytic continuation of 
I jy{uj) p to negative frequency is an even function and (ii) z^(a;) 7^ for all 
nonzero frequencies. The first assumption is needed in order to extend the 
frequency integrals to the negative real axis, while the second one ensure that 
all the reservoir oscillators couple to the system. 

For a purely magnetizable medium Xe = 0, and we can make f/ = by 
choosing the Coulomb gauge (V ■ A = 0). Let us define M = aX and write 
H as 



Later we prove that H in Eq. OlSI) is in fact the same H which appears in 
Maxwell equations. 

Now we go to reciprocal space and write all the fields in terms of their 
spatial Fourier transforms. For example electric field is written as 



(tilde distinguishes the field in real and reciprocal spaces). The total La- 
grangian in reciprocal space is 




(14) 



H = KoB - M. 



(15) 




(16) 



L 



f 




(17) 
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Since F (k,t) = F(— k,t) is satisfied for any arbitrary real field we should 
restrict the integration to half of the reciprocal space, which is indicated by 
a prim over the integral in f|T71) . The Lagrangian densities in the reciprocal 
space are obtained as 

Cem = eo{E^-c^B^), (18) 

C^at = P^^ - PUJI^^ (19) 

Cres= / duipY^-pU^Yl), (20) 
JO 

Cint = ak X A* ■ X - / duu{uj)X ■ Y* + c.c. (21) 

Jo 

The Coulomb gauge in reciprocal space can be written as k ■ A(k, t) = 0. 

Using the definition of longitudinal part of matter field, that is k x x" (k, t) = 
0, and the coulomb gauge condition, it is easy to prove that k x A(k, t) ■ 

x" (k, t) = 0. This shows that there is no interaction between the longitudinal 
part of the matter field and EM field. Therefor without losing any generality 
we can consider matter and reservoir fields as transverse fields. 

We introduce unit polarization vectors eA(k), A = 1,2, which are orthog- 
onal to k and to each other, and decompose the transverse fields along them 
to get 

A(k,t)= iA(k,t)e,(k) (22) 

A=l,2 

with similar expressions for the other fields. £ can be used to obtain the 
conjugate components of the fields 

dC ~ 

- eoEx = ^ = eoA^, (23) 
dAl 

Px = ^ = pX,, (24) 
dC ~ 

Q.x = ^— = pY^x - v{u)Y^x- (25) 

dy:x 
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Using the Lagrangian (fT7|) and the expressions for the conjugate variables 
we obtain the Hamiltonian as 

-f^ = y" rf^k(7l!em + Ti-mat + Ti-int), (26) 

where 

Km = £o(E)' + £ocVA'. (27) 

and 



p Jo p 

+ rrf^(!Mx*.QJ + c.c., (28) 
Jo p 

is the matter part, including the interaction between the magnetization and 
the reservoir. lUq = lUq + duj '"^^2 is the renormalized frequency of the 
polarization field and 

Hint = -a j d?k[\iyi k* ■y. + c.c], (29) 

is the interaction between the EM field and the magnetization. The fields are 
quantized in a standard fashion by demanding ETCR between the variables 
and their conjugates. For EM field components we have 

[iA(k,t),ll,(k',t)] = Z!^<J,;,,5(k-k'), (30) 
and for the matter fields 

[iA(k,t),P;'(k',t] = ^/i<5A,A'<5(k-k'), (31) 

[y^A (k, t) , g:,,, (k', t)]=ih 6xy 6{k - k') 6{u - u') , (32) 
with all other equal-time commutators being zero. 
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Using the inverse Fourier transform, the obtained Hamiltonian in (1261) 
could be written in coordinate space as 



H = / AlfggM + '"'^ " i^<^- "'^ V X A(r. t) . M(r. t)] 



2 

+ Hjnat- (33) 

The commutation relation between E(r, t) and A(r, t) from, fl30|) . is obtained 
as 

[A,(r,t),E,,(r',t)] =z-53!,,(r-r'), (34) 

where S^y {r—r') is the transverse delta function [14]. Using the commutation 
relations ( 134|) and the Hamiltonian ( l33l) . we find the Heisenberg equations 
for A(r,t) and E(r, t) as 

E(M) = %^, (35) 

and 

£0^^^ = KoV X V X A(r, t) - V X M(r, t). (36) 

Eg. 0351) is one of the Maxwell equations and Eg. (136 p is also the Maxwell 
equation in a magnetizable medium if we accept relation f|T5|) and interpret 
M(r,t) as the magnetization. This shows that the assumed interaction be- 
tween EM field and matter field lead to the proper Maxwell equations in a 
magnetic medium. 

In the following instead of solving the Heisenberg equation for EM field 
we write the EM field in terms of eigenoperators of the Hamiltonian (!26l) . 

To facilitate the calculations, we introduce a set of three annihilation 
operators as 



ax{k,t) = ^l^^[kcA^{k,t) - iEx{k,t)], (37) 



bx{k,t) = ^^jCo,Xx{k,t) + ^A(k,t)], (38) 
bx{k,u,t) = yi£[-z^X,,(k,t) + -^QUKt)]. (39) 
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From the ETCR for the fields, dSOD-dSD, we obtain the ETCR for above 
operators as 

[a^{k,t),a{{k',t)] = 6xx'S{k-k'), (40) 
[6,(k,t),6i(k',t)] = 5AA'5(k-k'), (41) 

[h{k,uj,t),b{,ik',Lo',t)] = 6xx'S{k - k')6{uj - to'). (42) 

We emphasize that, in contrast to the previous ETCR between conjugate 
fields, which where correct only in half space, Eqs. fl40l) - P2l) are valid in the 
hole reciprocal space. By inverting fl37|) - fl39|) to find the field operators, and 
inserting these fields into the Hamiltonian fl2B]) . we obtain after integration 

^ern = f d^k a{ (k) dx{k), (43) 

A=l,2 



Hmat = fd^^Y. {^^0 b{{k)bx (k) + r ujhuj h{ (k, u) hx{k, uo 

•' A=l,2 

+ -J^ duViu)[b{i-k) + bxmb{i-^,u;) + b,{k,u)]} 



(44) 



H^nt = fd'k Y: «\/^(«A(k) + d{{-k))J-^{by{k) + bU-k))e,y, 
-I X A'=l 2 V ^^0 V ^P^O 

(45) 

where eAA' is the antisymmetric symbol, cuk = c|k| and V{uj) = [^^] 



The polarization and reservoir parts of the Hamiltonian Hmat (1451) can be 
diagonalized by using the Fano technique [TB] to get a dressed matter field. 
The diagonalized expression for Hmat is (the detail of this technique can be 
found in [5] and we only give the results) 



Hmat= r duo I d^'k J2 huBi{k,u)Bx{k, 

^0 J A=l,2 



Lu), (46) 
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where i?^(k, a;) and i?(k, cj) are creation and annihilation operators of the 
dressed matter field respectively, which satisfy the usual ETCR 

[Bx{k,uj),B'^{k',uj')] = 6xx'S{k-k')6{uj-uj'). (47) 

They can be expressed in terms of the initial creation and and annihilation 
operators as 

B{k,uj) =ao{u)bik)+po{uj)b\k)+ / dujaiuj,uj')b{k,uj)+P{uj,uj')b\k,uj), 

Jo 

(48) 

and all the coefficient ao{u), f3o{uj), ai{uj,Lj') and /5i(u;,a;') can be obtained 
in terms of microscopic parameters. 

Using the commutators of b with B and B'^ together with fj^Hj) . it is easy 
to show that 

bx{k)= dLo[al{uj)Bi{k,u)-[io{uj)B,^{kuj)]. (49) 

We call this model the damped magnetization model. In this model the 
explicit form of ao(w) and (3q{uj) in terms of the microscopic parameter such 
as matter field density, p, or the coupling between matter and reservoir, 
is not important and we can only accept that such a function exist. 

In terms of the new set of operators the interaction part of the Hamilto- 
nian can be written as 



//■oo h X 

(fk dujY.-K{k)g{uj)Bi{k,u){ax{k) + a\{-k))exx' (50) 

•^0 A,A' 2 

where A(fc) = and g{uj) = («*(.;) - P*{uj)). 

Now we can follow two different methods to obtain the time dependence 
of EM field. One is writing the Heisenberg equation and solving it by using 
the Laplace transformation. The second one is using the Fano technique and 
finding the diagonalized H and writing EM field in terms of new operators. 

Here we following the second method and write the Hamiltonian as 

H= r duj I d^k V nwC{{k,uj)Cx{k,Lo), (51) 

■^0 A=l,2 

where 
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Cx(k,Lj) = ao{k,uj)ax(k) + (]o{k,u)a{ik) 

+ / du' [a{k,u,uj')By{k,uj,uj') + P{k,uJ,uJ')Bl,{k,uJ')]exx'■ 
■^^ A'=l,2 

(52) 

If we follow the algebra of reference [5] the coefficients in the relation 
above, fl52l). can be written as 



- ujlz{ujy 



ao{k,uj) = { — - — )— — (53) 



7, X ^ cj-g^k ^ Viuj,k) 



«(.,.,.') ^ Si. - + (f )(^^^)(^), (55) 



where V{uj, k) = A{k)g{uj) and z{uj) = 1 - /_+^ ^^'^^ and e 0+.^ 
Following the method used in fHOj) to derive the matter field operator b 
in terms of the dressed matter operators B and B'^ , we invert 0521) to write 
the photon annihilation operators a and dressed matter operator B in terms 
of the new operators C and C'^ as 

ax{k) = r duj[a*{k,uj)Cx{k,uj) - Po{k,uj)Cl{k,uj)], (57) 







Bx{k,uj)= duj'[a*{k,uj',u)Cx(k,uj) - f]{k,uj',uj)CUk,uj)]. (58) 

JO 

From Eq. fl37|) . A(r, t) is given by 



A(r,t) = -1^ / rf^k Y: J^[dxik,t)e^^-^ + H.c.]exik), (59) 
(27r)2 ^ ^^^ 2 V 2eoa;k 
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and we can use this expression and relation to obtain 



fu^l r,^_ /(-) 



(27r)2 i ;^tl^2 V 2£o -'0 u;2 - a;^ + 
X [C'A(k,a;)e'('^-^-'^*) + //.c.]eA(k) (60) 



where /(o^) is defined as f{uj) = — and Xm(w) is defined by 

2 J-oo a;' — a; — 2 J-oo — w 2 



(61) 



Later we will show that Xm('^) in Eq. (!6T|) is the magnetic susceptibility. It 
can be seen from fICT]) that the obtained magnetic susceptibility satisfies the 
Kramers— Kronig relations. In addition f{uj) can be written as 



|/(a;)r = ^^!^^. (62) 

TT 



Let us write fl60l) as 



A(M) ^ .(^) /A E r,.^^^n^p^ 

X [C'A(k,cj)e*(^-'-^*) -/7.c.]eA(k). (63) 

Relation fl63l) only depends on magnetic susceptibility which is a macro- 
scopic quantity. This result is exactly the same as the result of macroscopic 
approach in a magnetizable medium [T^ . 

Now we show that Xm('^) is the Fourier transform of the magnetizability 
of the medium. For this purpose we write 



X(r,t) = J d^k Y: [h{k)e^^- + H.c.]e,{k), (64) 



A=l,2 



and using relations ( l49l) and ( l58l) we easily obtain 



- / Pj f + CO 



X' 



X Cy{k,u')ex'x}- (65) 
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Using fl60p . Eq. fl65p can be written in real space as 



Jo 



'OO 



M(r,t) 



du{[KQxU^)V X A(r,a;) + M;v(r,c<;)]e-*'^* + i^.c.}. (66) 



Comparing fl^B]) with constitutive equation in ([HD, the interpretation of 
x(a;) as a magnetic susceptibihty is confirmed. In addition, the magnetization 
noise operator in ([8]) is obtained as 



M;v(r,u;) = [ d'k J2 pnc'eoImXm{i^)Cyi\^,uj)e,^^,e'' \ (67) 



Considering ([S7D, and the ETCR between M(r,cu) and m\Y,uj') 



[Mmx{j,uj),M\^{v',uj')] = 2tichoImXm{uj)5xx'6{v - v')6{uj - J). (68) 



Relation (l68l) satisfies dissipation— fluctuation theorem and is the same as ex- 
pression of noise operator in macroscopic method [12] . Therefor the supposed 
Lagrangian is equivalent to macroscopic approach. 

4 Extension of the Huttner-Barnett model to 
a magnetodielectric medium 

In pervious section we showed that the HB model can be extended to a 
magnetizable matter by changing the interaction between EM and the matter 
fields. The extension of this model to a magnetodielectric matter can be done 
by considering the following Lagrangian density 



A'=l,2 



becomes 




(69) 



where Cem is defined in ( ITT]) , and 



C 



Imat 




(70) 



P^2 pu;^ 2 
2 ^ 2 ^ 



(71) 
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are the matter parts of Lagrangian, which are modeled by two distinct har- 
monic oscillator fields Xi and X2 with the same frequency Uq. In the following 
we call Xi and X2 the polarization and magnetization fields respectively. It 
should be noted that taking the same frequency for both polarization and 
magnetization fields does not affect the result. The Lagrangians describing 
the reservoir are defined by 

C,res= r^cuf^YL-^YLV (72) 



Cres = I dcO (^|YL - ^Ytj , (73) 

and _ 

poo 

Clint = -aiA ■ Xi - dt^i^i(a;)Xi.Yi^, (74) 

Jo 

C2int = Oi2V X A ■ X2 - / da;z/2(u;)X2.Y2a;, (75) 

Jo 

are the interaction parts. The coupling functions i'i{uj) and 1^2 (i^)? satisfy the 
same assumptions we assumed for v{uj) in f|T^ . 

The displacement and magnetic fields are defined respectively by 

D(r,t) =£oE(r,t) + P,(r,t) (76) 

and 

H(r,t) = KoB(r,t)-M(r,t), (77) 

where we have defined P(r, t) = — aiXi(r,t) and M(r,t) = +a2X2(r,t). 

Again we choose the Coulomb gauge, but since the medium is polariz- 
able the scaler potential does not vanish. Since tJ does not appear in the 
Lagrangian, U is not a proper dynamical variable and should be written in 
terms of the proper dynamical variables A, X, and Y^. This can be done 
by going to reciprocal space and writing the Lagrangian as 

= / C^^k(-£em + ^ -^imai + -^ires + -^iint)- (78) 
i=l,2 

The Lagrangian densities in this space are obtained as 

Cem = eo& - KoB\ (79) 
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tl^at = - PO^I^I. (80) 

= / ^^a;(pY,, - pa;^YL) (81) 

with i = 1, 2, and 

= [-«iA* • i - /" (ia;zy(a;)X • "^^l + c.c, (82) 

>C2int = X A* • X2 - / (ia;i/(a;)X2 • Y^^ + c.c. (83) 

JO 

i * 

Using the Euler-Lagrange equation for U we find 

(/(k,i)=,i!i(lLXM). (84) 

Now if we decompose the longitudinal and transverse parts of the fields, then 
the total lagrangian can be written as the sum of two independent transverse 
and longitudinal parts as 

L = lII+L^, (85) 



where 



and 



= /' S^L^,^ + i^Lat + C,i,s + J , (86) 

■' i=l,2 



r 2 



>CL = eo(A -c^B^), (87) 

'*-imat — P^i ~ P^O ) 



= /^°° cIu;(pyL - ^Yi,2) , (89) 



= -ai A -k^- r du;i^^{u)X^ ■ Y^* + c.c, (90) 
Jo 

= a2k X A • X^^ - r du;i^,{iv)Xi ■ + c.c. (91) 
Jo 
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As in the previous section, we put away the longitudinal part of magneti- 
zation field. So the longitudinal part of the Lagrangian consists of only the 
polarization and EM field. Using fl84|) . fl79l) and flHOj) the longitudinal part of 
the Lagrangian in terms of the polarization field can be written as 



'ci^k4, (92) 



where 



£}[ = pXf - pulxf + f duj{pYll - pa;2Y£) 

Jo 

C+oo _ II _ II _ II _ I 

- / d^i;i(a;)(x!*-YL + x!*-YL). (93) 

Jo 



The is the longitudinal frequency and is defined hy ujl = \Jujq + uj^ where 

2 

ojI = The link between the transverse and longitudinal parts is given by 
the total electric field, which is written as 



~ A'l ~ II 

E(k, t) = -A(k, t) + — X"(k, t). (94) 

Using (p^ and the definition of the displacement field D given by (176|) . we 
recover the fact that the displacement vector is a purely transverse field. 

In this work, we are mainly interested in the transverse fields and shall 
only present the detailed quantization of the transverse part of the La- 
grangian. In what follows we restrict ourselves to transverse fields and omit 
the superscript ±. We use the same unit polarization vectors eA(k), A = 1, 2 
an find the the conjugate variables from C as 

dC ~ 

- soEx = = eoAx, (95) 
dAl 

Pix = = P^ix - (96) 

dC ~ 

P2X = = pX2x, (97) 

dX^x 
dC ~ 

Q±ujX = — — = pYi^x - v±{uj)Y±^x- (98) 
f)Y* 
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For the particular type of the couphng between hght and polarization 
field (?), the conjugate of A is the transverse electric field — eoE. A canonical 
transformation leading to a E ■ X type of coupling, gives the displacement 
field — D as the conjugate of A. Naturally, these two possibilities lead to 
the same results. We choose here the first possibility in order to keep as 
close as possible to the classical theory, where E is usually considered as the 
fundamental variable. 

Following pervious section we can obtain the Hamiltonian from the La- 
grangian (!86l) and the conjugate variables (p5!) - (!98l) . as 

H= f d^T^em + E ^i^at + (99) 
i=l,2 

where 

Hem = eoiEf + eoulA\ (100) 
is the electromagnetic energy density and (Dk is defined by u;^ = c^(A;^ + k"^) 
with kr = — = the Hamiltonian densities 



~2± ~± 2 

p Jo p 

+ rdu{^xt*-Qt) + c.c., (101) 
Jo p 

are the energy densities of the matter fields, where V^^uj), {i = 1,2), are 
defined in f H5|) . and 

7i:i,„ = — (A*. Pi + C.C.), (102) 
Pi 

Q/r, 

n2^n = -{kx A* ■X2 + C.C.), (103) 

P2 

are the interaction energies between the EM field and the matter fields. 

As usual, we demand ETCR between the variables and their conjugates 
and introduce the annihilation operators as 



dxik,t) = J^^{kcA,{k,t)-tEl{k,t)), (104) 
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kx{k,t) = L^{uoX,,{k,t) + -P*,{k,t)), (105) 

&i,(k, ,u;,t) = ,l^{-zujY,^^{k,t) + -^Ql^^{k,t)), (106) 

which satisfy the standard bosonic commutation relations. 

The Hamiltonian fl99|) can be written in terms of the annihilation and 
creation operators defined in Eqs. (11071) — (11091) . Since the structure of the 
matter fields in the Hamiltonian (llOip are the same as Hmat in (EHD, so the 
Hamiltonian of the matter can be diagonalized using the same method. The 
total Hamiltonian in terms of eigenoperators of the matter and reservoir fields 
can be written as 



H = y"(i2k{;k::;kai(k)aA(k) + ^°°(icj;kj^S|;,(k,a;)4A(k,a;) 

+ -Ai(fc)/ du;{g,{u;)Bl{k,u;)[a,{k) + d{{-k)] + H.c.} 
2 Jo 

- -A2{k) duY,{g2{uj)Bly{k,u)[ax{k) + ai{-k)]exx' + H.c.}}, 
2 Jo ^ 

(107) 

where the annihilation operators of the polarization and magnetization fields 



are 



b,x{k,t)= du[ao{uj)BUu,k,t)- (3o{uj)B,x{uj,k,t)], (108) 
Jo 

where A,{k) ^ , A^ik) = g^{uj) = ^(aoi(^) + /5oi(^)), 

(72 (cj) = (aQ2(^) ~ Po2i^)) the k integration has been extended to full 
reciprocal space. 

In Eq. (ll07p . EM field is coupled with two distinct reservoirs and is differ- 
ent from the usual Huttner model which only contains one reservoir, so we 
represent here the details of the diagonalization process of the Hamiltonian 

(unzD. 
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The diagonalization of H can be achieved by introducing the operators 
(^(k, u) as 

C'A(k,a;) = ao{k,uj)ax{k) + Po{k,uj)a{{k) 

+ / du'[ai{k,uj,uj')Bi{k,u') + (3i{k,uj,uj')Bl{k,uj')] 
Jo 

+ / du'y][a2{k,u,uj')B2X'{K^') + P2{k,u,u')Bly{k,uj')]ex.x', 

Jo y 

(109) 

where the coefficients are chosen such that the operators C(k, cj) satisfy the 
eigenoperator equation 

[C{k,uj),H] = hujC{k,u). (110) 

This equation, together with the expansion of the Hamiltonian (11071) and 
the definition of C(k,u) in (11091) . lead to the following linear equations be- 
tween the coefficients 



aQ{k,u)u! = ao{k,uj)ujk 
1 

+ 5^-/ duj'[a,{k,u,u')V,{k,u')-P*{k,u,uj')V,{k,u')], 

fllll 



(5o{k,uj)uj = -(3o{k,uj)uJi, 
2 Jo 



\ roc 

+ Eo/ du'[a,ik,uj,uj')V,ik,u')-P,{k,u,uj')V,ik,io% 



;ii2) 



a±{k, u, uj')lo = -[ao{k, u) — (3o{k, u!)]V*{k, u') + q;i(A;, oj, uj')uj' , (113) 
2 



(5^{k,u, Lo')uo = ^[ao(A;, u) - (3o{k, uj)]V*{k, to') - (3,{k, u;')u', (114) 



where Vi{k,u!) = Ai{k)gi{ 
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These set of equations can be easily solved to obtain f3o{k, u), a±{k, u, u') 
and (3i{k, u, u') in terms of ao(k, uj). Subtracting flll2p from fillip we obtain 

(3o{k,uj) = ^ ^^ ao{k,u). (115) 

We now replace ffTTSD for /3o(k,a;) in ffTT3|) and ffTTD . and find 

a,{k, u, Lo') = [P{^-^) + ydk, io)]V*{k, u;')^^ao{uj), (116) 

i3,ik,u,uj') = [^^]V,ik,u')^^aoiKu;), (117) 

CO + UJ' U + U]f_ 

where P means the Cauchy principal value. The relation between functions 
yi{k,uj) and y2{k,u>) can be obtained by substituting the expressions for 
ai{k,u>,u') and (3i{k,uj,uj') in flll6p and flll7p into fillip . Using the defini- 
tions of V^{k,u) and V2{k,u in fllOTp it is easy to show that they are odd 
functions of frequency uj. We use this fact to extend the integral in the nega- 
tive frequency region and obtain the relation between functions yi{k,uj) and 
y2{k,uj) as 

V,\k,u)yr{k,u) + Vi{k,u;)y2{k,u;) = "^^-^ + ]-P dJ^i^^ 

CJk 2 J-oo UJ — UJ 



2 J-oo U' — UJ 

(118) 

In order to calculate ao(A;,a;), we impose the standard commutation re- 
lation on (^(k, u) 

[Cik, uj), ^^(k', uj)] = 6{uj - u')6{k - k'). (119) 

Using the expression for C(k,uj) given by (11090 and the set of equations 
ao{k,u) flllSp — flll7p . we can find ao^k^u) (up to a phase factor) in terms of 
the yi{k,uj) and y2{k,u). By taking a suitable phase factor and doing some 
routine but tedious calculations we find the following expression for ao{k,u) 



cDfc \y,{k,uj)+trryV,^{k,uj) + {y2ik,u) + t7ryVi{k,u)^"- 



(120) 
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From the Eqs. fllU9p and fillip — ( TTT^ . we can obtain two independent sets 
of operators, C and C' which satisfy the following commutation relation 

[C{k,uj),C'\k,uj)] = Q. (121) 

For obtaining these operators, we first choose yi{k,u) and y2{k,u) such that 
they satisfy flllSp and then using fll20p . we find C and . For obtaining C 
and C'\ we should choose y[ and y'2 such that they satisfy (11 180 and also the 
derived operators C' and C"''" from them should satisfy (11210 . So, for defining 
y[ and 2/2 there are two equations and accordingly they can be determined 
uniquely, operator C. 

Since in Eq. (lll8p . yi{k,uj) and y2{k,uj) can not be determined uniquely, 
so we have a freedom in determining operators C, C' . But we do not lose any 
generality by taking a special solution since these operators are all equivalent 
up to a Bogoliubov transformation. 

To facilitate the calculations, we choose yi{k,uj) and y2{k^uo) such that 
they satisfy in (11180 and the following relation 

Vi{k,u)[yi{k,uj)-i'K] = +V2{k,u)[y2{k,u)-ni]. (122) 

Therefore yi{k,u) is obtained as 



yiik,^) = T7271 V—T7T, ^TTT, H ^+0/ duJ 



Vi{k,iL!) + Vi{k,uj)V2{k,uj) 2 J- 00 u' — u 

+ 0/ duj— i7cV2{k,uj){ . -I}- 

2 u' — u V2[k,uj) 

(123) 

Using (I120p and (I123p . we find ao{k,u>) as 

where Zi{k,u) = ^ /q°° -^t^r^- The other set of operators, C"(k, cj), can be 
obtained from 

Vi{k,uj)[y[{k,uj) — itt] = —V2{k,uj)[y2{k,u>) — irc]. (125) 
For aQ^k^u) we have 
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Eqs. (11241) and fll26p can be used to obtain C and C in terms of a, a^, Bi 
and S|. From ^Mj, and ([m]) we can write the Hamiltonian fllU7p as 

H = jd^kj^ hw[C\k,uj)C{k,u;) + C'^{k,uj)C'{k,uj)]. (127) 

Using the commutation relation (I119P and the commutation relation between 
a and a^, we can invert the Eq.( lll9p to write a and in terms of C, C\ C' 
and C"''" as 

a(k) = / dLu{a*Jk,Lo)C{k,uj) - Po{k,uj)C\k,uj) 
Jo 

+ a'*{k,uj)C'{k,u)- P'o{k,uj)C'^{k,uj)}. (128) 

Before writing the EM field in terms of C and , for later simplification 
we use a Bogoliubov transformation as 

^ C{k,u;) + &{k,u) 

Ke{k,u) = , (129) 

^ C{k,co)-&{k,u) 

K^{k,u) = -j= . (130) 

Using fl59l) and (I128p . A can be obtained in terms of eigenoperators of the 
Hamiltonian, as 



A(r,t) = ' ^[d^kHduil- -^^H^^(k,a;) 



(27r)fV2£o^ Jo ''ujl{l-xM)-^\'^ + Xe{uj)) 

^k/2Hi^m(k,a;) i^*(k.r-c.t) _Hc} 



(131) 

where 

XeM-^pd^'-%^ = ^prrf-'^ + ^H/iHP, (132) 

2 J-oo uj — uj — le 2 JO a; — a; 2 
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and 

Xm^ = o / ^ = / 7+o«7r /2Hh 133 

2 J-oD 00 — uj' — te 2 JO cij — cij' 2 

and = ^cDoC^fc^^i^ and ^(a;) = ^Mal^. 

Now the process of calculating M in the previous section can be repeated 
for P and M (P and M are defined in (1761) and fl77|) respectively). We find 

P(r,f)= / MhXeHE(r,a;) + P^(r,cu)]e-*"* + i/.c.}, (134) 

JO 

and 

POD 

M(r,t)=/ c?a;{[«:oXm(^)V X A(r,cu) + Mjv(r,cu)]e-*"* + //.c.}, (135) 

JO 

where PAr(r,a;) and MAr(r,u;) are 

pNxir,u) = I d'\^phe^ImXeKex{K^)e''''\ (136) 



Mnx{j,uj)= fd^k ^2heocHmXmK^x'{K i^)e'''"'exy. (137) 

•' A'=l,2 

By comparing f ll34p . f ll35p and (IH]), (E]) we find that Xe and Xm are electric 
and magnetic susceptibilities. As in the preview section, the commutation 

relation between P7v(r,u;) and p]v(r,i^) and M7v(r,u;) and M^(r,ijj) can 
be calculated. The results are compatible with the dissipation-fluctuation 
theorem and coincide with macroscopic results. 

Using relation ([IS]), (IIS2D and Eq.([ISSD, we find 



1 X /■ ,3, /""^ , rr uj\/ImXeiuj)Ke(k,u) 



A(r,t) = -^( ) d'k du{[ 



7r%o J Jo a;^(l -Xm(w)) -a;2(l + xe(w)) 



V/mx^(u;)i^^(k, g;) 
a;^(l-x™(^))-^2(l + Xe(^))^ ■ 



(138) 



The relations (11381) . fll32p and (I133p . and the commutation relations be- 
tween the noise operators are exactly the same as the results obtained from 
the macroscopic method [121 US]- So, these two methods are equivalent. 
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5 conclusion 

The Huttner-Barnett model has been extended to a magnetodielectric medium, 
The results obtained in the present model are equivalent with those obtained 
in the phenomenological models. The exphcit form of the noise operators 
have been obtained. Based on the results obtained here, the Lagrangian in- 
troduced in the present work can be used as a microscopic model for canonical 
quantization of the electromagnetic field in a magnetodielectric medium. 
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